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Abstract 



We present the analytical calculation of the contribution of order 
m^a^nj to the absorptive part of the vacuum polarization function of 
vector currents. This term constitutes an important gauge-invariant 
part of the full 0{rriia^) correction to the total cross-section of e + e~ 
annihilation into hadrons. The results are compared to predictions 
following from various optimization schemes. 
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1 Introduction 



In QCD the correlator of two currents is a central object from which important 
physical consequences can be deduced (for a detailed review see, e.g. jTj). In par- 
ticular, important physical observables like the cross-section of e + e~ annihilation 
into hadrons and the decay rate of the Z boson are related to the vector and 
axial- vector current correlators. Furthermore, total decay rates of CP even or CP 
odd Higgs bosons can be obtained by considering the scalar and pseudo-scalar 
current densities, respectively. 

From the theoretical viewpoint the two-point correlators are especially suited 
for evaluations in the framework of perturbative QCD (pQCD) |2j. Indeed, due 
to the simple kinematics (only one external momentum) even multiloop calcu- 
lations can be analytically performed. Non-perturbative contributions can be ef- 
fectively controlled through the the operator product expansion [3J H] . As a con- 
sequence, the results for practically all physically interesting correlators (vector, 
axial-vector, scalar and pseudo-scalar) are available up to order a 2 s taking into 
account the full quark mass dependence jSl El Ej- 

In many important cases (with the Z decay rate a prominent example) the ex- 
ternal momentum is much larger than the masses of the (active) quarks involved. 
This justifies to neglect these masses in a first approximation which significantly 
simplifies the calculation. As a result, in massless QCD the vector and scalar cor- 
relators are analytically known to a\ jHJEllIIIj- The residual quarks mass effects 
can be taken into account via the expansion in quark masses. At present this has 
been done for the quadratic and quartic terms to the same order [TTJ [T2J EH . 

During the past years, in particular through the analysis of Z decays at LEP 
and of r decays, an enormous reduction of the experimental uncertainty (down 
to 0(lCr 3 )) and in R(s) has been achieved with the perspective of a further 
reduction by a factor of four at a future linear collider [Ti] . 

Inclusion of the 0(a^) corrections jS] is mandatory already now. Quark mass 
effects as well as corrections specific to the axial current [T^l EH HZ] must be 
included for the case of Z-decays. The remaining theoretical uncertainty from 
uncalculated higher orders is at present comparable to the experimental one pQ. 
Thus, the full calculation of the next contributions, those of 0(aj), to R(s) is on 
agenda. 

In massless approximation R(s) is conveniently written as 




^2 Q) 3 f 1 + a s + of (1.98571 - 0.115295 n f ) + 



a^(-6.63694 - 1.20013 n f - 0.00517836 n)) + a 4 s r t 




(1) 
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where a s = a s (fi 2 — s)/tt and the standard MS renormalization prescription JS| 
is understood. The a 4 term can be further decomposed as a polynomial in rif, 
namely 

V,4 V,4 . V,4 . V,4 2 i V,4 3 

r o = r o,o + r o,i n/ + r o,2 n f + r / 3 n f . (2) 
The term of order a^n^ (and, in fact, all terms of order a s (a s nf) n ) have been 
obtained earlier by summing the renormalon chains ^H]. These are technically 
very simple to compute 2 , numerically small (r ' 3 = 0.02152) and not directly 
sensitive to the nonabelian character of QCD, in contrast to the terms of order 
ajnj. The first nontrivial contribution — the subleading 0(n 2 j) piece in (j2J) - 
have been only recently analytically computed in [20j with the result 3 

<2 4 = -0.7974. (3) 

Taken by itself eq. (JBJ) is not of much use for the phenomenology for obvious 
reasons. However, it gives a strong extra support to the well-known prediction 
of the full 0(af) contribution of the Kataev and Starshenko|2"T] (for a detailed 
discussion see [22]). The prediction is often used to estimate the theoretical uncer- 
tainty in R(s) due to higher order not yet computed perturbative contributions. 

Another point of concern are the quark mass effects at order a; 4 . Fortunately, 
simple estimates (see e.g. [TT] ) show that, say, the quadratic mass effects are 
completely under control 4 at the scale of the Z-boson mass since terms of order 
0(al) are available and small. This holds true for the vector correlator with 
the leading term proportional a s m1 as well as for the axial vector correlator, 
where the the leading in a s term is present already at the Born level and thus 
proportional to m 2 . Quarks mass effect are getting progressively more important 
at lower energies. From the conceptual viewpoint it would be important to test 
the evolution of the strong coupling as predicted by the beta function and the 
standard QCD matching procedure through a determination of a s from essentially 
the same observable, however, at lower energy. The region from several GeV above 
charm threshold (corresponding to the maximal energy of BEPC around 5.0 GeV) 
to just below the B meson threshold at around 10.5 GeV corresponding to the 
"off resonance" measurements of CESR or B-meson factory is particularly suited 
for this purpose. As a consequence of the favorable error propagation, 

a s\ 1V1 Z> 

the accuracy in the measurement (compared to 91 GeV) may decrease by factor of 

about 3 or even 4 at 10 and 5.6 GeV respectively, to achieve comparable precision 

in Aqcd- 

2 At least for a fixed number of loops. 

3 For brevity we display it in the numerical form. 

4 Provided one uses the properly chosen parameterization in terms of the running quark mass. 
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Technically speaking, the evaluation of the leading in n/ term of order 0(m 2 ajn 3 j) 
is again rather simple while the subleading term of order 0(m 2 a 4 n 2 ) presents 
a problem comparable to that in the massless limit. In this work we describe 
the corresponding calculation and its results. We limit the discussion to the 
vector current correlator, which is relevant for hadron production through the 
electromagnetic current. We also compare our results with predictions following 
from various optimization schemes. Our basic conclusion is once again that the 
FAC/PMS optimization predictions are remarkably close to reality. They provides 
a quantitative argument in favour of the corresponding full prediction. 



2 Generalities 

To fix notation we start considering two-point correlator of vector currents and 
the corresponding vacuum polarization function (j* = Q^^Q; Q is a quark field 
with mass m, all other n/ — 1 quarks are assumed to be massless) 

IW?) = iJdxe^(0\T[ 3 ;(x)f u (0)}\0} = (-g fll/ q 2 + q II q u )U(q 2 ). (4) 
The physical observable R(s) is related to IT(g 2 ) by 

R(s) = 127r3n(g 2 + ze). (5) 
It is convenient to decompose R(s) into the massless and the quadratic terms 

R(s) = 3 {tf + ^} + --- = 3jE^ 

Here we have set the normalization scale fi 2 — s; a s — a s (s)/ir and dots stand for 
corrections proportional to m 4 and higher powers of the quark mass. 

For the calculation of r^' 4 we had to deal with divergent parts of five-loop 
diagrams and finite parts of the four-loop ones (see, e.g the corresponding dis- 
cussion in 0I2S])- In fact, as was first discovered in Ref. [TT] . the quadratic mass 
correction r^' 4 can be obtained with the help of renormalization group methods 
exclusively from the four- loop function n 2 defined as (a s = a s (/i 2 ) ) 

n = fn + ^n 2 ) , n 2 = ]T a\k v 2 >\ 

16tt 2 V -g 2 J ^ 

Indeed, due to an extra factor of m 2 the function II2 obeys an unsubtracted 
dispersion relation. As a result the following homogeneous RG equation holds: 

U^- 2 + lm{a s )m^- + P(a s )a s ±) Tl 2 = 0. (6) 
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Equivalently, since IT 2 is independent of m, 

d ( d \ 

— n 2 = f -2 lm (a s ) - (3(a s )a s —\ IT 2 , (7) 

2 

where L = In Once the rhs of eq. (jJJ) is known, one can find the function 
n 2 up to a constant contribution which has no effect on R(s). Thus, to compute, 
say, the 0{— aj|) correction to R(s) one needs to know 

• The three- loop QCD (3 function function (starting at 0(a 2 )) and the quark 
mass anomalous dimension (starting at 0{a s ))\ available from each 
up to order a^. 

• The very polarization operator IT 2 at three loops (that is of order a 2 s ): was 
computed also long ago in 26J. 



Altogether, this implies that the 0(—af) term in R(s) could have been been 
computed as early as in 1986. But in reality this was done only by 5 years later 
[TJ. 

From the above discussion it is clear which ingredients are necessary to com- 
pute quadratic mass terms in R(s) at C(af): 

1. The four -loop QCD (3 function function and quark mass anomalous dimen- 
sion: available from [23 EB1 EH] • 

2. The polarization operator Il 2 at /our loops, including its constant part. 



3 Results 



Using the new technique described in j3THl3*H l2Ti] and the parallel version of FORM 
[321 we have computed the leading and subleading (in nf) contributions of 
order to the polarization function n 2 (0(a 2 ) results are known from |26j ) 



n 2 
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y a s + a s 



95 



n f + 



18923 784 „ 4180 , 
■ Cs + -z=- Cs 



54 



+ at 



5161 8 



' Cs 



1458 27 
20 , 28880 , 

+ — C4 - —T7T C5 



27 

62893 424 
+ 



27 



162 



27 



Cs 



4150 
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3 243 
or, numerically, 

n 2 = -8-21.333 a s +a 2 s (10.56^-224.80)+^ (-3.896^+274.37^+^), (9) 
where we have set a s = a s (— g 2 )/vr, jj? = —q 2 . 
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The straightforward use of eqs. ()5I7|) yields (a s = a s (s)/7r,/i 2 = s) 



a* 
12— 

7T 



Numerically 
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+ 2 

S3 



™/ + r 2 i o r V + r yi j 



rX = Yla s + a 2 (-4.3333 ra/ + 126.5) + a 3 (l. 2182 nj- 104.167 n/ + 1032. 
+ (-0.20345 n) + 49.0839 nj + r^ 4 n/ + r^' 4 ) • 



(10) 

14) 
(11) 



4 Discussion 

Let us contrast eq. (JBJ) with predictions based on the optimization methods FAC 
(Fastest Apparent Convergence) and PMS (Principle of Minimal Sensitivity) 34, 

EH! EI] • As is well known in many cases FAC and PMS predictions are quite 
close |SZ|- To compare with an independent prediction we will also use so-called 
NNA (Naive NonAbelianization) approach of [38J. 

In Table 1 we list the FAC and PMS predictions for the function n 2 (see <— 
|40j ) . From the three entries corresponding to n/ = 3,4 and 5 one easily restores 
the FAC/PMS prediction for the nj dependence of the a\ term in LT2 

fc^ 3 (FAC) = -4.2^ + 277^- 2886, (12) 
fcJ^PMS) = -4.2 nj + 282 n/ -2916. (13) 

The comparison with © shows very good agreement with the calculated terms of 
order nj and rA. As a natural next step we will now use the available exact infor- 
mation about the coefficients fc 2j ' 2 and k 2 \ to find k 2 ' from FAC/PMS prediction 
for a selected value of nf = 3. The choice seems to be natural as in many cases 
FAC/PMS predictions made for n/ = 3 are in better agreement with the exact 
results (see. e.g. j^U E3 1221)- We obtain (in fact, very close values are produced 
with rif = 4 and 5) 

fcg(FAC, nf = 3) = -2880, fcg(PMS, nf = 3) = -2897. (14) 

Now we turn to NNA. The corresponding prescription is to simply make the 
substitution n/ — > — 3/2/3 = nj — 33/2 in the leading term of the n/ expansion. 
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The result is 

fc^' 3 (NNA) = -1060.69 + 128.56871/ - 3.896 n). (15) 

This prediction is in obvious conflict with the result of the direct calculation of 
the term linear in rif. 

Table 1: Estimates for the coefficient k 2 ' based on FAC and PMS optimization 
schemes. 



Method 


71/ = 3 


77/ = 4 


77/ = 5 


FAC 


-2092 


-1844 


-1604 


PMS 


-2109 


-1856 


-1612 



5 Conclusions 

We have presented the analytical calculation of contributions of order m^ajn'j to 
the vacuum polarization function of vector currents. Predictions of FAC, PMS- 
and NNA optimization methods are tested against the exact results. Good quan- 
titative agreement is found for the first two cases. NNA seemingly predicts only 
sign and order of magnitude of analytical results. 
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ported by the DFG-Forschergruppe "Quantenfeldtheorie, Computeralgebra und 
Monte- Carlo- Simulation" (contract FOR 264/2-1), by INTAS (grant 00-00313), 
by RFBR (grant 01-02-16171), by Volkswagen Foundation and by the European 
Union under contract HPRN-CT-2000-00149. 
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